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Abstract: We present a new program (iHixs) which computes the inclusive Higgs boson 
cross-section at hadron colhders. It incorporates QCD corrections through NNLO, real and 
virtual electroweak corrections, mixed QCD-electroweak corrections, quark-mass effects 
rS ' through NLO in QCD, and finite width effects for the Higgs boson and heavy quarks. iHixs 

C^ \ can be used to obtain the most precise cross-section values in fixed order perturbation 

theory in the Standard Model. In addition, it allows for a consistent evaluation of the 
cross-section in modified Higgs boson sectors with anomalous Yukawa and electroweak 
interactions as required in extensions of the Standard Model. iHixs is interfaced with 
the LHAPDF library and can be used with all available NNLO sets of parton distribution 
functions. 
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1. Introduction 

A major objective of experimental high energy physics is the discovery of the Higgs boson. 
Production cross-sections at hadron accelerator experiments are expected to be small. 
However, they are sufficiently significant in order for the Standard Model (SM) Higgs boson 
to be discovered soon in a variety of signatures. A remarkable progress in this direction 
has been made after the direct searches of LEP [1]. The TEVATRON [2] experiments 
have already demonstrated sensitivity to Higgs boson cross-sections of magnitudes as in 
the Standard Model. ATLAS and CMS have published limits on the inclusive cross-section 
using first LHC data [3,4] accumulated in 2010. These searches are expected to yield tight 
constraints in their forthcoming updates with data from the 2011 runs. 

Interactions in the Higgs sector of the Standard Model are only indirectly constrained 
experimentally. Theory-wise, a modified Higgs sector is a prerequisite for extensions of 
the Standard Model with a more appealing UV completion. When computing Higgs bo- 
son production cross-sections, it is necessary to allow for the possibility of Higgs boson 
interactions with modified couplings from their Standard-Model values. 

At a hadron collider, the inclusive Higgs boson cross-section is most sensitive to the 
values of the quark Yukawa couplings. The main hadroproduction mechanism is gluon 
fusion which receives sizable contributions from top and bottom-quark loops. In Standard 
Model extensions, the mass of elementary particles may not originate entirely from the 
Higgs mechanism and, in addition, not yet discovered quarks may be postulated which 
also contribute to the cross-section. Very enhanced or very suppressed Higgs boson cross- 
sections may be obtained, as it occurs for example in a Standard Model with a fourth 
generation [5-8] and composite Higgs boson models [9-11] correspondingly. 

Direct production from quark partons in hadrons (mainly the bottom and charm) 
may also become sizable in models with more than one Higgs doublet [12] or dynamically 
generated Yukawa couplings [13, 14]. In models with a perturbative Higgs sector, gluon 
fusion is likely less sensitive to modifications of the Higgs boson couplings to electroweak 
gauge bosons; these start contributing to the production of a single Higgs boson only at 
the two-loop level. The bulk of the generally small electroweak corrections is due to loops 
with light Standard Model quarks which have very well constrained electroweak couplings. 

A light Higgs boson is expected to have a small decay width in the Standard Model. 
The width could be different in extensions of the Standard Model. For example, the Higgs 
boson could have a significant invisible decay width in scenaria with hidden sectors [15-17]. 
The width of the Higgs boson is expected to be large also for scenaria where a heavy mass 
is not disallowed. 

The objective of this publication is to provide precise theory predictions for the produc- 
tion cross-section of a Higgs boson in the Standard Model and extensions which alter the 
interactions of the Higgs boson and quarks or electroweak gauge bosons. This publication 
is accompanied with a computer program, iHixs, for the precise calculation of inclusive 
Higgs boson cross-sections at the Tevatron and the LHC in such a general setup. 



2. Features of iHixs 

In the last two decades many theoretical studies lead to improved estimates of the Standard 
Model Higgs boson total cross-section at hadron colliders with a level of 10 — 20% precision. 
Using, extending and combining the available theoretical calculations has been proven to be 
a non-trivial task for both theorists and experimentalists. iHixs is an easy to use program 
which produces accurate predictions for the Higgs cross-section in the Standard Model and 
can be adapted readily for a large class of extensions of the Standard Model. 

For the numerical evaluation of loop and phase-space integrals we have used a new 
Fortran package [18] of harmonic polylogarithms with complex arguments. For one-loop 
box and triangle master integrals with different internal masses we have used the library 
OneLOop of Ref. [19,20] which allows for complex masses. We have tested our implementa- 
tion of all other one-loop master integrals against both OneLOop and the numerical program 
of Ref. [21]. 

iHixs computes the inclusive Higgs boson cross-section through next-to-next-to-leading- 
order (NNLO) in perturbative QCD. A large source of uncertainty for Higgs boson cross- 
sections at hadron colliders is the precision in the determination of the parton densities. 
It is therefore important to compare the effect of diverse existing determinations of parton 
densities on the Higgs cross-section, as well as future sets which will incorporate refined 
measurements and theory. iHixs allows these studies effortlessly. It is interfaced through 
the LHAPDF library [22] with all available parton distribution functions with a consis- 
tent evolution at NNLO [23-25]. Other sets of the library can be employed by simple 
modifications of the code. 

iHixs allows the study of the Higgs boson invariant mass distribution for a finite 
width of the Higgs boson and compute the cross-section sampling over a Breit-Wigner 
distribution. This assumes that matrix-elements for the production of a Higgs boson decay 
factorize in matrix-elements for the production of a Higgs boson times matrix-elements for 
the decay. A grid of values for the branching rations and the decay width of the Higgs boson 
in the mass range of interest is a necessary ingredient for iHixs. For the Standard Model 
(or models with very similar Higgs decay rates), iHixs uses a grid of values for the width 
and branching ratios that we produced with the program HDECAY version 3 . 532 [26]. To 
study models where the width and branching ratios differ significantly from the SM, the 
user has to provide his own grid, in the form of a simple text file. 

iHixs includes perturbative contributions for two production processes: 

- gluon fusion 

- bottom-quark fusion 

2.1 Components of the gluon fusion cross-section in iHixs 
The cross-section for the gluon fusion process in iHixs comprises: 

1. leading order and next-to- leading order (NLO) QCD effects with exact quark-mass 
dependence. The number of quarks and their Yukawa couplings are arbitrary. 



The required two-loop amplitude has been first computed in Ref. [27,28] where it was 
presented in the form of an integral representation ^ . This was later expressed in terms 
of harmonic polylogarithms with the method of series expansion and resummation in 
Ref [31]. Independent analytic evaluations were performed in Refs [32,33]. 

The real radiation matrix-elements have been computed in Ref. [34, 35] and recom- 
puted for the purposes of several other publications including this one. Numerical 
implementations of the NLO QCD cross-section with full quark-mass effects in the 
Standard Model were made in Ref. [28] and in Refs [36,37]. 

2. NNLO QCD corrections, using heavy quark effective theory (HQET). 

The NNLO Wilson coefficient for an arbitrary number of heavy quarks and Yukawa 
couplings has been computed in Ref. [11]. In the special case of Standard-Model 
Yukawa couplings this is equivalent to the Wilson coefficient of Ref. [38] , while in the 
case of integrating out only a single heavy quark it yields the Wilson coefficient in 
Refs [39,40]. 

The NNLO phase-space integrated matrix-elements in HQET have been first eval- 
uated in their threshold limit [41,42]. The complete NNLO correction has been 
computed in Refs [43-45]. 

3. two-loop electroweak corrections at leading order in a. 

These include the full Standard Model contributions to the amplitude as computed 
in Ref. [46,47] ^. For a light Higgs boson, a dominant contribution to the two- loop 
amplitude is due to loops with light quarks [48]. We allow for a common re-scale 
factor of the HWW and HZZ couplings. This should be a sufficient parameter for 
models which have a custodial symmetry protection in order to comply with stringent 
constraints from electroweak precision tests ^ . 

4. one-loop electroweak corrections for the real radiation processes qq — )■ gh and qg — t- 
qh. 

This amplitude has been first computed in Ref [49]. In this paper, we perform an 
independent calculation and present analytic formulae in terms of a basis of finite 
master integrals adding also contributions with massive quarks in the loop. Our 
results agree with the limit of zero Higgs boson mass of Ref. [49]^. Subprocesses with 
bottom quarks in the initial state have also been studied in [80] at a differential level. 



^ First calculations of the two- loop amplitude in the infinite top-quark mass limit were performed in 
Refs [29,30] 

•^ We thank the authors of Refs [46, 47] for kindly providing a text file with the numerical values of the 
two-loop amplitude. 

^Custodial symmetry protects from large corrections to the T parameter; see, for example, the constraints 
on custodial symmetry breaking operators of the strongly interacting Higgs boson effective Lagrangian 
in [10]. 

*We were unable to compare with the analytic expressions for the electroweak amplitude of Ref. [49] 
due to the lack of an exact definition of the contributing "finite parts" from the divergent master integrals 
which were chosen as a basis. 



5. mixed QCD and electroweak contributions with light quarks. 

This contribution can be estimated by means of an effective field theory and the 
required Wilson coefficient computed in Ref. [50] . 

6. arbitrary Wilson coefficient for the HtT{G^yG^^) operator. 

New physics at scales higher than the electroweak scale may introduce modifications 
to the gluon fusion cross-section which cannot be accounted for by modified Yukawa 
interactions and rescaling electroweak corrections. In such situations, iHixs allows 
to introduce corrections to the Wilson coefficient of the HQET theory. Given the 
non-discovery of new states with Tevatron and first LHC data, it is reasonable to 
anticipate that the new energy frontier is distant enough from the electroweak scale 
(higher than the top-quark mass) in order for an effective theory approach to be 
adequate. 

2.2 Components of the bottom-quark fusion cross-section in iHixs 

While the dominant production mode of the Higgs boson in the Standard Model is gluon 
fusion via a top-quark loop, a considerable correction of about 5% arises from bottom-quark 
loops. Another production channel based on bottom-quark fusion is gg — >• Hhb [51-53]. In 
this channel the bottom-quarks are predominantly produced collinear to the gluons. The 
cross-section at fixed order in perturbation theory suffers from large logarithmic terms. 
However, these can be resummed into the bottom parton density, leading to another 
(mostly) single Higgs-boson production mechanism bb ^ H [54-56]. We have therefore 
included this channel into iHixs. 

The cross-section for bottom-quark fusion process in iHixs comprises the NNLO cross- 
section calculation of Harlander and Kilgore in Ref. [57]. This cross-section is important 
when the Yukawa couplings to bottom quarks are enhanced, as it may happen in models 
with more than one Higgs doublet. We have implemented the analytic formulae of Ref. [57] 
and computed the scale-dependent terms of the cross-section separately. We have checked 
that our numerical code agrees with the publicly available program of Ref. [58]. 

We included the bottom-fusion process in iHixs for the purposes of facilitating the 
simultaneous study of enhanced bottom Yukawa couplings in the production of a Higgs 
boson from bottom-quarks and in gluon fusion via bottom-quark loops. We note that 
the iHixs program can be adapted easily in order to compute the cross-section for Higgs 
production via the fusion of lighter quarks, such as the charm-quark, if necessary [14]. 

3. Higgs boson interactions 

We consider a Higgs boson with interactions described by the Feynman rules: 






SM 




SM 



The triple- vertex Higgs-quark-quark is the product of an arbitrary, flavor dependent factor 
Yf and the analogous Feynnian rule in the Standard Model. iHixs permits an arbitrary 
number of quark flavors Nf in order to accommodate extensions of the Standard Model 
with novel quarks. 

The Standard Model Feynman rules for the H — W — W and H — Z — Z vertices 
are rescaled by a global factor Xewk- We did not find it necessary to introduce a separate 
re-scaling factor for the W and Z boson vertices. The ratio of the coefficients of the 
corresponding operators is fixed by the custodial symmetry and very tightly constrained 
by electroweak precision tests [10]. 

4. Hadronic and partonic cross sections 

We consider the production of a Higgs boson, H^pn), at a hadron collider, which decays 
to some combination of final state particles that will be collectively denoted by {.fffinai}: 

hadroni(PL) + hadron2(P2) -> {H(p//) + X, (other processes)} -> {-f^final} + X (4.1) 

The hadronic cross-section is given by the factorization theorem as, 



ijSpartons 

where 



J2 j dxidx2fi{xi,iif)fj{x2,iif)al^^^jj^^^^^^^{s,iif) (4.2) 



full _ \^ I ^^.rl^^fJ^. ,..\f.(^^ ,,.^;:^/"" 

^{J^final}+^ 



XIX2S, s={Pi + P2y. (4.3) 



The indices i,j run over the flavours of initial state partons. The functions fi{x,fif) are 
parton distribution functions in the MS— factorization scheme and Hf is the factorization 
scale. 

Singling out typically dominant contributions from resonant diagrams as p^ — )■ m'jj, 
we cast the partonic cross-section in the form 

'■full ^ _. '■signal— bkg _. ^bkg (A A\ 

The first term on the right hand side corresponds to the square of the resonant Feynman- 
diagrams, the second term corresponds to the interference of resonant and non-resonant 
diagrams, and the last term is the square of diagrams without a resonant Higgs propagator. 
We shall refer to the first term as the "signal cross-section". It can be written as: 



Q^h{Q) ^ij^nis, Q% fif)BTH^{Hf,^,,}{Q) 



Qa,Qb define the experimentally accessible range for the invariant mass of the particle 
system originating from the decay of the intermediate Higgs boson, mn is identified with 
the Higgs-boson mass. Away from the resonance region, Q^ ~ mj^ the above equation is 
adequate for the signal cross-section. iHixs is dedicated to the evaluation of the "signal 
cross-section" . A full description of the cross-section for a Higgs final state can be obtained 



by adding to the "signal cross-section" the remaining two contributions of Eq. 4.4: a theo- 
retical or experimental estimate of the background cross-section and a theoretical estimate 
of the signal-background interference. 

In the resonant limit, Q^ — )■ w-^, the signal cross-section is dominant and it becomes 
infinite at any fixed order in perturbation theory, for Q^ = m^ exactly. A resummation 
of resonant contributions at all orders is necessary in order to render the propagator finite 
in this limit. We remark that a resummation of partial perturbative corrections from all 
perturbative orders into the propagator of an unstable particle is a delicate theoretical 
issue [59,60]. Historically, it has been treated with various prescriptions in the literature 
with varied success (see, for example, references in [61]). To a first approximation, the 
signal cross-section becomes: 

F/f (Q) is the decay-width of a Higgs-boson at rest with mass Q. In the zero Higgs boson 
width limit it reduces to the product of the partonic production cross-section dij^n for an 
on-shell Higgs boson times the branching ratio for its decay ^XH^ffj^^_^^\. Performing the 
transformation 

— ^ = 1 -I- otan(7ry), o = (4.7) 

m|^ rriH 

we obtain an equivalent integral with a better numerical convergence, 

cvb QTiiiG) 
^iJMHfin.i}+xis,f^f)= J dy—^^^—^dij^H{s,Q'^,^if)Bi„^{H^^^^}{Q). (4.8) 

The integration boundaries 7/^5 are computed from Eq. 4.7. 

A light Higgs boson, as predicted in the Standard Model, has a rather small 5 and it 
is often sufficient to take the 5 = limit of the zero width approximation (ZWA). Existing 
experimental studies at hadron colliders [2-4] have always reported limits on the Higgs 
boson cross-section comparing with expectations in this approximation. However, recent 
years have witnessed the alarming trend of using this approximation in situations where it 
may be severely insufficient^ . We therefore find it useful to dedicate a part of our numerical 
studies to Higgs bosons with a non-negligible width. 

The width grows for heavier Higgs boson masses due to decays into electroweak gauge 
bosons. In extensions of the Standard Model, this feature may be more pronounced as new 
decay modes may be available. For such situations, the zero width approximation is poor. 
An integration over the Breit-Wigner distribution of Eq. 4.8 is a more accurate estimate of 
the signal cross-section. We also note that Eq. 4.8 convolutes with the branching ratio. For 
Higgs masses close to thresholds branching ratios are steeply changing; a naive estimate of 
the signal cross-section in the zero width approximation could also be unsafe. 

iHixs allows the possibility for the calculation of the Higgs signal cross-section taking 
into account finite Higgs width effects by performing the integral of Eq. 4.8. Notice that 



^for example, recent limits on the cross-section for a heavy Higgs boson [3] 



the resummed expression of Eqs 4.6,4.8 is in good agreement with Eq. 4.5 away from 
the resonant region only when the partial width, T h (Q)^^ H^iHa^^AiQ) 7 is computed at 
a variable Higgs-boson virtuality Q. Higgs decay rates are rather sensitive to the Higgs 
boson virtuality due to the many decay mechanisms which become available at diverse mass 
values. In order to use iHixs with an arbitrary BSM model, the user needs to provide 
a data file with the width and branching ratios of the Higgs boson as a function of the 
virtuality of the Higgs boson. 

For the Standard Model, or models with similar enough width and branching ratios, 
we have generated a grid of values using the program HDECAY of Ref [26] ^. We note 
that for large Higgs boson widths the description of the Higgs line-shape may require fur- 
ther improvements, both in the resonance region where a more sophisticated resummation 
framework could be employed [59] and, especially, for virtualities far from the Higgs boson 
mass where a dedicated estimation of the signal-background interference cross-section is 
needed. We believe that iHixs can provide a flexible enough platform for such modifica- 
tions if need arises (e.g. with experimental evidence of a Higgs boson and the associated 
heavy particles which render the theory consistent with electroweak precision tests and 
unitarity bounds). 

In the Standard-Model, it has been observed that significant cancelations due to 
interference of resonant and non-resonant diagrams take place at high invariant masses 
(Refs [62-65]). The magnitude of the "signal-background" cross-section is very important 
and cannot be neglected. iHixs takes into account only diagrams with an s-channel Higgs 
boson propagator. The line-shape away from the resonance is therefore poorly described. 
To improve upon this, we have implemented a prescription based on the resummation of 
VV —7- VV scattering amplitudes with the dominant contributions from both resonant and 
non-resonant Feynman diagrams [66] at the high energy regime. The contributions which 
unitarize the scattering amplitude for vector-boson scattering can be approximated with 
the amplitude for Goldstone boson scattering at high energies. Conveniently, the amplitude 
in this regime can be described in terms of an effective Higgs propagator. Ref [66] performs 
a Dyson re-summation of the tree-level Goldstone boson scattering amplitude leading to an 
"improved s-channel approximation" . In this framework, the Higgs propagator is modified 
according to the prescription: 






.mt 



m 



(4.9) 



H s - mjj + iT H{mjj] 



rriH 



This prescription interpolates smoothly between two limits which are well described 
either by resummation or by fixed-order perturbation theory: the resonant region Q '^ mn 
and the high energy limit Q S> rriH- We do not envisage Eq. 4.9 as the final step towards 
a precise description of the line-shape for the Higgs boson. However, it is a very useful 
diagnostic tool in order to assess how important the signal-background interference could 
be for a heavy Higgs boson. According to this prescription, the hadronic cross-section is 

^For studies with iHixs where the HDECAY tabulated Higgs boson width and branching ratios are used 
we request that Ref [26] is also cited 



computed as, 



x/seym(Q,m-H), (4.10) 



with 



^ _ m 

I seym "^ 



4 ii-<r\\g2 



^'•'^-^ <l^,^r <"'' 



™h/ '^h 



We emphasize once again that when Eq. 4.10 is used (Seymour option in iHixs) some 
signal-background interference effects which are dominant at very high invariant masses 
are taken into account. In contrast, Eq. 4.8 (default option in iHixs) computes purely 
the signal cross-section (only the square of resonant Feynman diagrams). 

We believe that the Higgs boson line-shape will enjoy many future theoretical studies 
with improved resummation methods for resonant diagrams and matching to fixed-order 
perturbation theory away from the resonance region. In the course of these developments, 
new prescriptions which are not yet implemented in iHixs shall emerge. We have made 
efforts to be able to include such improvements readily in iHixs. For example, a method 
introduced recently in Ref. [67] requires that all virtual amplitudes must be computed 
with a complex Higgs virtuality Q^. This requirement is effortless to achieve in iHixs 
where we evaluate all QCD one and two-loop amplitudes using CHAPLIN [18] for harmonic 
poly logarithms with complex arguments. We defer to the future a conceptual appraisal and 
a numerical comparison of theoretically appealing approaches beyond the prescriptions of 
Eq. 4.8 and Eq. 4.10. 

We now define the dimensionless ratios 

r=^ f = ^ z^^ (4.12) 

s s xia;2 

The hadronic Higgs signal cross-section can be cast in the form 



^-^ [y^ qth{Q) p .^. , ,^ 

JJSpartons J'" 

} ^ / — dzCij{xi,X2,lJ-f) 



X 

where 

T 
Xl =X, X2 = , Cij = [Xifi{xi)] [X2fj{x2)] . (4.14) 

xz 
In iHixs we have implemented the partonic cross-sections for two Higgs boson production 
channels: 

(i) bottom-quark fusion as in Ref. [57] 
(ii) gluon fusion. 

Amplitudes for the two processes do not interfere and the two cross-sections can be com- 
puted independently. 



We have faithfully adopted the definitions and analytic expressions for the partonic 
quark-fusion cross-sections of Ref. [57], where they have been computed through NNLO 
in the strong coupling expansion. We shall therefore not discuss this process any further, 
except for the presentation of numerical results with iHixs. 

The perturbative evaluation of the gluon fusion cross-section has been the topic of nu- 
merous publications. Nevertheless, for the sake of clarity, we find it important to elaborate 
on our implementation of perturbative corrections for this process. 

5. The gluon fusion process through NLO QCD 

We cast the partonic cross-sections in the form 



z 288^2^ V TT 

where [ir is the renormalization scale for the strong coupling. 

(v) 
The partonic cross sections zrv^- at NLO and higher orders consist of virtual and real 

emission parts that are separately infrared divergent. We expose the infrared singularities 

of the real radiation matrix-elements by the method of plus-distribution subtractions. We 

then add the virtual part, the collinear counter term which has been generated from parton 

distribution factorization in the MS— factorization scheme and the renormalization counter 

terms in the MS— renormalization scheme. The resulting cross-section is finite and can be 

written as the sum of three distinct terms: a term proportional to (5(1 — z) that corresponds 

to all contributions from leading order kinematics (virtual part plus 5-proportional terms 

from the integrated soft-collinear pieces), a 'regular' term that corresponds to higher order 



kinematics, and plus-distribution pieces, proportional to 
We therefore write 



1-2 



for various f(z). 



nif {Q'^,z,tif,fj,r) = ^ n'^^i^{Q^,z,fif,fir) (5.2) 

ke{5,+,R} 
5.1 LO: gg ^ h 

At leading order only the gluon gluon subprocess contributes. One obtains 






Y^Y,t,^A{t,) (5.3) 



with 



Tg - ^2 ^^-4^ 

and niq being the mass of the heavy quark in the gluon fusion loop. The quantity Tq^A{Tq) 
has the simple limits 

lim Tq-AiTq) = 1 lim r,-^(r,) = (5.5) 

Tq— S-OO / Tg— >U Z 



The full analytic expression is given by, 



1 fl + r l^ 



(5.6) 



with 



(yr^^+1) 



(5.7) 



Note that the sum runs over all quarks in the model, and Yq is as defined in the Feynman 
rules of Section 3. 

5.2 NLO:gg-^h + g 

The (5-part of the NLO correction to the gluon gluon subprocess can be written as 



e-li^P 



2/3olog(|]|)+vr2 



+ K 



B^V, 



(5.8) 



where 



B^Y.^,r,-A{T, 



(5.9) 



is the LO coefficient, and 



o q21 

7,{t,) = Y,-Mfl{T,) = Y,^ 



(5.10) 



r(l), 



where MrJ^{Tq) can be found in eq.7.4 of ref. [32] and Q^ can be found in eq. 26-30 of 
ref. [33]. The plus distribution part is 






-6.og|| 



1-Z 



+ 12 



log(l - z) 
\- z 



(5.11) 



Finally, the regular part of the NLO gluon gluon correction is 

2 



,(1) 



"'99;R 



(1-.)A(1-A)\2 ^ 



+|5|^U 



Vgg{z)\0% 



{\-zf 



log 



\-z 



X-z^z^fm^ 



/"f 



61og(§^)p..(^)K5-12) 



with Pgg{z) the gluon splitting kernel 



Pggi^) = - + z{l- z)-2 



(5.13) 



The form factors ^i^q„H{Tg) can be found in the Appendix B.2. 
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5.3 NLO: qq^h + X 

The qq initial state starts contributing to the total cross section at order a^ in QCD. There 
are also non-negligible mixed QCD-electroweak corrections, at order a^aewk which we shall 
discuss later. The pure QCD corrections lead to the following coefficients: 

(5.14) 
(5.15) 



n 



(1) 



^% = n 


(1) _o 


32 (1 - zf 

27 z 


/ ^^qTqAqqgH{z) 

q 



qq;R 27 

with Aqqgniz) given in eq. B.5. 

5.4 NLO: qg^h + X 

The gluon-quark initial state also contributes at order a^, and it receives similar mixed 
QCD-electroweak corrections of order a^aewk- The pure QCD coefficient is: 



,(1) 



(1) 



%g;5 - ^qg 



0, 



(5.16) 



and 



n 



(1) 
q9;R 



\B\ 



-Jiz-Pgg{z)\Og^^^_^^^ 



tJ'F 



+ 



dX 



with 



(1-A)_ 



yx 



Yl ^q'^q^qqgH ivx 



Pgg{z) log ( ^ 

'l + (l-z)2A 



—z 




(l-z)(l- 


A) 


Cf 1 + (1 


-zf 



Pgqiz) 
and AgggH^z) given in eq. B.5. 

5.5 NLO: mixed QCD-EW corrections to qq ^ H + g and qg ^ h + g 

The mixed QCD-electroweak contributions^ have the following structure: 



(5.17) 

(5.18) 
(5.19) 



n 



{0)ewk 
qq;R 



1 



Kwh"^- "Y "Y ^\ TqAqqgH{z)cx,q ' 

q Xe{Wi,Zi,H} [ 

2 2 

17* ^ ^"*13 I 7-1* / \ '^23 



Sl2 



Sl2 



(5.20) 



where Aggguiz) given in eq. B.5. The sum over X runs over all W and Z-like bosons in the 
model, as well as over the Higgs boson, while the sum over q runs over all heavy quarks, 
6, t, . . ., as before. The coupling cx,q contains Kronecker delta symbols which select specific 
initial state quarks, depending on X, as explained below. 



''Note that these corrections are of order a^, so we denote them by n, 



(O)e'wk 
qq 
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• Z loops: For the initial state quarks q G {u, d, c, s, 6} we have the following couplings 
\z = 2 



CZ,q = {5qu + 5qc){Vz^u + ^zj + (V + ^qs + 5qb){Vz,d + ^Z, 

I'z,^ = ^(tt - ^sm fc*^) , az,u 



d> 



cos9,„ 2 3 ' 2cos^. 



w 



fl'"' / I ■ 2 /) \ ffw /'COIN 

^^z,d = Z-(-o + Q '^^ ^"' ' "'Z,d = t; s- (5-21) 

cos t^^ 2 3 2 cos Uw 

The form factor is identical for all initial state quarks and yields 

-^i" = -inlAewk {s3i, S23, si2,Q,mz) 

Jg" = -mlAewk is23, S31, si2,Q,mz) (5.22) 

with Aewkis,t,u,mh,mz) given in the Appendix in eq. B.7 

W loops: The couplings are given by 

Xw = 2, cw,q = ivw,q + aw,q)^ '^w,q = —7=-, aw,q = —7^, (5.23) 

J^r = -'miAewk{s3i,S23,Si2,Q^,'m^)^SqL-'ml,A]^^l,{s3i,S23,Sl2)Sqb 

L 

TY = -m^yle«,fc (s23,S31,Sl2,Q^,m^) ^^gL -m^^^^(s23,S3i,si2)M5-24) 

L 

where L sums over all light quark states u,d,c,s. Note that here we have summed 
over the internal light quark flavors. This yields a CKM coefficient of 

Y^ \Vijf ^1 for i = l,2 
i=i,2 

Measurements show this to be true to about 1 in 10000, so it is a good enough 
approximation to make. 

For q G {b}, we have to take the internal quark to be a top. The couplings are 
unchanged if we use that IV33P ~ 1, which is also a good approximation to make. 
The form factor A^^ {s, t, u) is given in the Appendix in eq. B.9. 

Higgs in the loop: In the standard model this gives a non negligible contribution for 
q = b, other quarks may be considered as well if their Yukawas are enhanced. The 
couplings are given by 

Xh = 3, CH,q = ^qb{v%^q + ojjg), VH,q = ^ ^ , aH,q = 0, (5.25) 



and the form factors are 



-T^f = -mlAH{s3i,S23,Si2) 

^2 = -mlAff (s23, S31, S12) . (5.26) 
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6. Beyond the NLO QCD 



The leading order and next to leading oder cross-section are known exactly. Beyond that 
we can only take the limit of heavy electroweak gauge bosons and top-quarks. Bottom 
quark contributions are also unknown beyond NLO. 

In the effective theory approximation, the Higgs gluon interaction is described by an 
operator of the form 



-Ceff 



1 






where the Wilson coefficient has a perturbative expansion 



(6.1) 






vr 



vr 



+ ... 



(6.2) 



The cross-section calculated in the effective theory is, 



C^eff 



Co + C,^iM + c/"^(^^ 



IT 



IT 



+ ... 



Vo H Vi + — V2 + 

vr V vr / 



(6.3) 



Expanding in as we obtain 



with 



0"cff 









"eS ^ ■ ■ ■ 



(6.4) 



O-eff = |Co|%0, 

a^^ = |Co|%i + 2Re(C7oCi)r/o, 



a 



(2) 
eff 



|Co|%2 + 2Re (CoCi) r/i + ( |Cip + 2Re (C0C2) ) r?o. 



(6.5) 
(6.6) 
(6.7) 



The integrated cross-sections rji have been computed through NNLO in Refs [43-45]. 

Resorting to an effective theory calculation is necessary only for important QCD and 
electroweak corrections which cannot be evaluated in the full theory. In the previous 
sections, we have listed results in the exact theory for the LO and NLO QCD perturbative 
expansion as well as one-loop electroweak corrections. In addition, two-loop electroweak 
corrections are also known exactly [47]. We will keep these corrections with their full mass 
dependence as in the exact theory calculations and use the effective theory approach for 
contributions at higher orders in the strong and electroweak couplings, namely for the 
NNLO correction in QCD and mixed QCD and electroweak corrections. 

We match the effective theory and full theory perturbative expansions as follows. Let 
us assume that we can compute the contributions to the cross-section exactly through 
some perturbative order for only some of the heavy particles which contribute to Higgs 
production amplitudes: 



_ (0) fas\ (1) /Os 

Impartial — (^partial '^ \ ^ ) partial '^ \ ^ 



partial ' ' ' ' 



(6.8) 
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In an effective tlieory approach, these contributions would factorize as in: 



C^partia^eff 



c, 



partial ^partial '^s (,/^j ^partial / 0!.s\fJ') 



+ C[ 



+ C' 



IT 



vr 



+ ... 



a. 



vo H — m + — V2 + 

vr V vr 



a.c 



We then write the cross-section as 



n=0 



^ = E V partial©!" < Norder) + 5a^^s 



(6.9) 



(6.10) 



where Norder is the last perturbative order that the partial contributions are known in the 
full theory and 



6a 



6a 



(0) 
cff 

(1) 
eff 



I Co I - 6(0 < Norder) 



a 



partial 





partial 



m, 



I Co I - 6(1 < Norder) Gq > r?i 

+ {2Re (CqCi) - 6(1 < Norder)2Re (^^P^rtiai ^partial ^ | 



6a 



(2) 
eff 



C, 



partial 



I Co I - 6(2 < Norder) C'J > r?2 

+ {2Re (CqCi) - 6(2 < Norder)2Re (^^P^rtiai ^partial ^ | 



m, 



m 



(6.11) 



(6.12) 



+ <^ |Cir- 6(2 < Norder) 



C'^ 



partial 



+ 



{2Re (C0C2) - 6(2 < Norder)2Re (^cP'-^t-'Cf ^""^) }) %. (6.13) 



6.1 The Standard Model Wilson coefficient with anomalous Yukawa and elec- 
troweak couplings 

In our theory, we can integrate out the top-quark, the new heavy quarks, and the elec- 
troweak gauge bosons W and Z. This yields a Wilson coefficient which is 



Co = Aqcd • 1 + Aewk • 1 
Ci = Aqcd • -J- + Aewk • ^ 
C2 = Aqcd • C2q + Aewk • C2w ■ 

The factor for the electroweak component is given by 



Aewk = A 



3a 



ewk 



IGtts?,, 



4 + 4 



* O ^^ A 

e^ J a^ 

4 3 "" 9 "' 



(6.14) 
(6.15) 
(6.16) 

(6.17) 



with Sw, Cyj the sine and the cosine of the Weinberg mixing angle. The factor for the QCD 
component is given by the sum of the anomalous Yukawa coupling re-scaling factors for all 
heavy quarks 

\qcd= Y. ^1- (6.18) 

gS heavy 
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The three-loop QCD Wilson coefficient C2q has been computed recently in Ref. [11] and can 
be read from Eq. (3.35) of the same reference. The four-loop electroweak Wilson coefficient 
is not yet known. 

C2w = unknown (6.19) 

One can attempt a rough estimation since it is conceivable that the perturbative series 
for the QCD component and the electroweak component follow a similar pattern: \C2w\ ~ 
\C2q\ ~ 10. In our studies, we vary C2w G [—30,30] as an estimate of the higher order 
mixed QCD and electroweak corrections. 

In iHixs we have implemented the exact contributions up to NLO, i.e. Norder = 1, 
for the heavy quarks and the electroweak gauge-bosons. These contributions should not 
be counted twice in the effective theory calculation. 

Cr'''^' = Aqcd • 1 (6.20) 

^partial ^ ^^^^ . H ^^^1) 

CP-t-i = Aqcd • C2q. (6.22) 

6.2 Improving on the effective theory approximation 

It has been observed that the effective theory works better for the K-factors rather than 
the absolute cross-section. At next-to-leading order, all real and virtual amplitudes in 
the soft or collinear limit have the same dependence on the masses of the heavy quarks 
as at leading order. It appears that the factorization of the cross-section in the infrared 
limit closely resembles the factorization of the cross-section in the limit of infinitely heavy 
massive particles. Finite quark-mass effects are important for "hard radiation" terms, 
but these are expected to have a typical perturbative expansion where an a^ suppression 
occurs from one order to the other. Top-quark mass effects have been studied with explicit 
calculations with operators of higher dimension in HQET demonstrating the validity of the 
approach at NNLO [70,71]. We can then improve on the effective theory approximation 
by making the replacement 

Aqcd^ Y. Yq^TqAiTq), (6.23) 

q G heavy 

and 

Aewk -^ KwkMggh^^v^K {mjj, M^, M|, . . .) (6.24) 

the two-loop electroweak amplitude for gg — )■ H, computed fully in [46,47]. In Fig. 21 of 
Ref. [46] we find the quantity 

2 



5ewk/100 = ^^^^+<3^^ 



—born 
"QCD 



-1 

top only 



1 I ■^ggh,EWK 



2TqA [Ttop) 



1. (6.25) 



We then substitute, 

3 

Aewk -^ Ae^fc-r^A (rtop) x Aq, (6.26) 
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with 



Ao 



1 + 



^EW 



K 



100 



1. 



(6.27) 



The two-loop electroweak corrections were kindly provided to us in a data file, electroweak . h, 
by the authors of Ref . [46] . 

7. Numerical results in gluon fusion 

In this section, we present numerical results for the Higgs boson cross-section via gluon 
fusion. We will first make a short discussion of the stability of the perturbative expansion 
and the scale variation uncertainty. Then we shall compare predictions from all available 
NNLO sets of parton distribution functions. We will proceed with a study of finite width 
effects for the heavy quarks in the gluon fusion loops. Finally we shall discuss the finite 
width effects on the Higgs boson total cross-section. To the best of our knowledge, there 
have been no published result for NNLO K-factors for the signal cross-section beyond the 
zero width approximation for the Higgs boson. 

7.1 Perturbative convergence and scale uncertainty 

The gluon fusion cross-section exhibits a rather slow convergence of the perturbative series 
in the strong coupling constant. iHixs computes the cross-section through NNLO in 
perturbative QCD. The perturbative behavior of the total higgs cross section with its scale 
uncertainty, when using MSTW08 parton densities, is shown in Fig. 1. for LHC at 7TeV 
collision energy. One notices that radiative corrections are sizable, where neither the NLO 
nor the NNLO corrections can be neglected. 
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Figure 1: Inclusive Higgs cross section at LO, NLO and NNLO, with scale uncertainty bands, 
calculated in the range n G [^o/2, 2/io] for MSTW PDFS. 
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Figure 2: Inclusive Higgs cross section at LO, NLO and NNLO, with scale uncertainty bands, 
calculated in the range /x e [mo/2, 2/io] for ABKM and GJR PDF sets. 

A major source of theoretical uncertainty in Higgs production via gluon fusion is due to 
the choice of the factorization and renormalization scales. We evaluate the scale uncertainty 
with a central scale of fip = fiR = fiQ = mn/^, and a variation in the range n G [Aio/2, 2^o]- 
We note that the NLO scale uncertainty band engulfs the NNLO band in almost the entire 
mass region depicted. We also notice that the magnitude of the perturbative corrections 
(K-factor) is larger for lower values of the Higgs boson mass. 

Similar behavior is observed when using the ABKM set, but not with the GJR set. 
Fig. 2, where the NLO and NNLO overlap is only partial. 

7.2 PDF comparison 

An enormous progress has been made in the last decade towards improving and estimating 
reliably the precision of parton densities. While uncertainties are generally small, the gluon 
fusion process requires the less constrained gluon density. The extraction of this quantity 
is an active field of research. It is very important to compare the effect of various gluon 
density determinations on the Higgs boson cross-section. We have enabled the possibility 
for such studies in iHixs. 

In tables 5,6,7 we present the inclusive cross section for pp — )• H + X for Higgs masses 
ranging from 110 to 300 GeV, using three different PDF sets available at NNLO, namely, 
MSTW2008 [24], ABKM09 [23] and GJR09 [25] . We have restricted our choice of pdf 
sets to the ones with NNLO DGLAP evolution, consistently with the iHixs computation 
of the partonic cross-sections through the same order. We believe that these iHixs results 
constitute the most precise predictions in fixed order perturbation theory for the Higgs 
boson cross-section. 

A comparison between the three different NNLO PDF sets and the corresponding 
PDF-|-as uncertainty bands are shown in Fig. 3. The uncertainties quoted include varia- 
tions of Os around the preferred value for every set, and are estimated according to the 
prescriptions of the PDF providers. A combination of cross-section values for a large set 
of pdf parameterizations is necessary, and these are efficiently computed in iHixs simulta- 
neously. 
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Figure 3: The Higgs production cross section at NNLO witli three different PDF sets and their 
68%CL FDF+fls uncertainty bands. 



The reason for the remarkably different predictions between the different PDF sets 
(that range from 10% in the low mass region to 30% in moderately high masses of around 
300GeV) is hard to trace. The bulk of it may be attributed to differences in the adopted 
values of the strong coupling constant as{mz)- The situation is only partially remedied if 
one chooses to consider the 90% CL (as opposed to the one-sigma, 68%CL) uncertainty 
bands provided by the MSTW collaboration. 

The comparison, in Fig. 4, shows the ABKM and MSTW uncertainty bands to marginally 
overlap. We finally note that preliminary results [68] with the updated ABMIO PDF fit, 
which includes hadron collider data, show that larger values for the total cross section are 
obtained, in comparison with ABKM09. 



7.3 Top quark width 

iHixs evaluates two and one-loop amplitudes in all kinematic regions, permiting a definition 
of mass and kinematic invariants in the full complex plane. This is a particularly useful 
feature when a resummation of finite width effects in threshold regions is necessary. For 
example, in a resummation framework using the complex mass scheme [60, 67] the masses 
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Figure 4: The Higgs production cross section at NNLO witli three different PDF sets and their 
PDF+Qfs uncertainty bands, using the MSTW 90% CL grids. 



of heavy quarks need to be evaluated according to the prescription, 



nig -^ niq {niq - iVg) 



(7.1) 



where Tg is the total decay width of the quark and niq its mass. 

Using iHixs, we have studied the finite width effects for quarks in fermion loops. We 
find that the top width is insignificant (at the level of less than one per mille) for a Higgs 
boson mass below the tt threshold. Around and above that threshold, its effect grows to 
the percent order as shown in Fig. 5. 

7.4 Finite Higgs boson width effects 

In Section 4 we discussed that there exist various approaches on how to treat the Higgs 
propagator when departing from the zero width approximation (ZWA). In this section we 
will present numerical results for the two different prescriptions described in Section 4: the 
default scheme (DEF) of Eq. 4.8 and the Seymour scheme (S) of Eq. 4.10. 

In Fig. 6 the inclusive Higgs cross section calculated within the ZWA and the two 
finite width schemes is shown, as a function of the Higgs mass. The width is calculated by 
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Figure 5: Relative difference Scrnw/crnw — " f • 100% of the cross section for the top quark with 
a real mass, cr*, and in the complex mass scheme with Ftop — 2 GeV. 

interpolating over a detailed grid^ constructed with HDECAY [26]. The cross sections^ are 
shown in table 1 . We note that the three calculations deviate widely for Higgs masses larger 
than 300GeV. The deviation between the ZWA and the finite width schemes is expected 
since for large Higgs masses the width of the Higgs boson is comparable to its mass. It is 
also evident that the finite width schemes deviate from each other in the high mass region, 
indicating a possibly large contribution due to signal and background interference which 
the Seymour scheme attempts to simulate. 

Within this context, it is interesting to notice that the invariant mass distribution 
of the Higgs boson, shown in Fig. 7, gets significantly distorted in the high mass region, 
where the Higgs width is large. The distortion is spectacularly stronger in the case of the 
Seymour scheme, as a consequence of the fact that the scheme tries to simulate the effects 
of signal-background interference off the resonant peak. These effects become increasingly 
important for high Higgs masses. 

In experimental searches for the Higgs boson where its invariant mass can be recon- 
structed from the momenta of the final state partons, as is the case for H — t- 77, or 



*The precision of the interpolation is always better than 2-10 ^ . 

^We use here the MSTW PDF set. Similar behavior is observed when using the other two NNLO PDF 



sets. 
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Figure 6: Comparison of tlie total cross section in tlie zero widtli approximation, o"^^"*, with a 
finite width in tlie default scheme, a^^^ and in the Seymour scheme, a^^^ . In the lower panel we 
show the relative error one makes when adopting the ZWA, defined as '^~zwa ■ 100% 
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Table 1: Total cross section for LHC at v^ = 7TeV with MSTW PDFs with a finite width in the 
two schemes, (jDEF,sey ^ ^^^ -^^ ^j-^g ^ero width approximation denoted by a^'^^. 



H — )■ ZZ, it is beneficial for the analysis to impose a kinematical cut on the total invariant 
mass of the Higgs decay products. The invariant mass of the Higgs boson is then con- 
strained in a window around the nominal Higgs boson mass, the size of which depends on 
the experimental resolution, see e.g. [72]. As a consequence, part of the signal is also cut. 
The signal cross section that survives such kinematical cuts on the Higgs boson in- 
variant mass can be calculated with iHixs and is shown in table 2, for the window sizes 
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Figure 7: The invariant mass distribution of the Higgs boson with mn = 200,400,600,800 GeV, 
in the defauh and the Seymour scheme. 

employed in [72]. We observe that the reduction in the expected signal rate can reach 
20% — 40% for window choices smaller than the Higgs width. A non-negligible reduction 
of a few per cent persists even when the invariant mass window is larger than the nominal 
Higgs width, due to contributions from the tail of the Breit-Wigner distribution. 
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Table 2: Total cross section, (jDEF,sey compared with the cross section in the invariant mass 
region tuh ± 5Q, denoted by a^EF.w ^^ ct^ey-w^ f^j. ljj^. at ^i = 7 TeV with MSTW PDFs . 

This effect can be estimated by parton shower Monte Carlo simulations which are the 
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main simulation tools in experimental collaborations. It is important that a realistic line- 
shape for a heavy Higgs boson is implemented in these simulations. Common practice in 
experimental studies is to evaluate distributions with a LO or NLO Monte Carlo program 
interfaced with parton showers and then rescale the distributions by inclusive K factors. 
Those K-factors, however, depend on the scheme adopted for the Higgs width, as well as 
on the size of the experimental window, if one exists in the analysis, as shown in table 3. 
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Table 3: NNLO K-factors with the width according to the Default scheme, K^^^q, in the presence 
of kinematical windows, K 



DEFy. 
NNLO ' 



within the Seymour scheme, Kf^^^Q^ and within the Seymour 
scheme in the presence of kinematical windows, Kj^j^£q, compared with the K-factors in the zero 
width approximation, -ft'^^VLO- 



When one departs from the zero width approximation, the branching ratios into the 
various final states also depend on the virtuality, as opposed to the nominal mass, of the 
Higgs boson. Assuming, at first, that all invariant masses are reconstructed experimentally, 
table 4 shows the difference between convoluting the branching ratio to a WW final state 
with the production cross section and the Breit-Wigner distribution, as in eq. 4.6, and 
just multiplying the total cross section with the branching ratio evaluated at the nominal 
Higgs mass value. We see that, in the low mass region, the relative deviations are very 
small, at the per mille level, thanks to the stability of the branching ratio in the region 
sampled by the Breit-Wigner. In the high mass region they can become large, especially 
when off-resonant effects are taken into account, as is the case for the Seymour scheme. 

An interesting theoretical question is whether the resummation recipe for the Higgs 
propagator employed in the vicinity of the Higgs resonance affects the tails of the invariant 
mass distribution, where we know that the correct propagator is the one appearing in 
eq. 4.4. To study this, we compare the invariant mass distribution with the distribution 
calculated in the region of the tails with the width of the Higgs boson set to zero in the 
denominator of the propagator. This distribution diverges at the peak, as expected, so 
we cut off a small region around the peak, of the order of r(?7i//)/2 or smaller, to assist 
convergence. The comparison is shown in fig. 8, for the default scheme and in fig. 9 
for the Seymour scheme. In both cases the tails of the distributions are described well, 
when \Q — m//| > T{mH) which indicates that the precise prescription for the propagator 
resummation in the peak region doesn't affect the tails of the invariant mass distribution. 
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Table 4: Cross section convoluted with the branching ratio to WW, compared with the product 

(Ti = (Jpp-fH X BRh-^ww, the relative deviation 6ai = ( — 1) • 100%, and the product 

with the production cross section in the ZWA, (72 = cr^Y^^ x BRh^ww^ with the relative deviation 



5CJ2 = (- 



pp^H^WW 



1) • 100%. All numbers are for LHC at ^i = 7TeV with MSTW PDFs 
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Figure 8: The invariant mass distribution of the Higgs boson with niH — 200,400,600,800 GeV, 
in the default scheme, compared to the distribution of the tails computed with the off-resonant 
propagator. 

7.5 Inclusive Higgs boson production in the presence of a fourth generation of 
quarks 

An extension of the Standard Model with an additional family of quark and leptons yields 
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Figure 9: The invariant mass distribution of the Higgs boson with itih — 200,400,600,800 GeV, 
in the Seymour scheme, compared to the distribution of the tails computed with the off-resonant 
propagator. 

a large increase to the gluon fusion process. It is therefore tested more easily at the 
TEVATRON and the LHC than the SM scenario . 

The most accurate computation of the total inclusive cross section for the production 
of the Higgs boson in a model with a fourth fermionic generation, for the LHC, has been 
presented recently in [69]. This calculation was performed in the ZWA. 

The value of the physical Higgs width in such a scenario depends on the details of the 
model, and is generally bigger than the Standard Model Higgs width. In this paper, we 
would like to assess the impact of the width, making the rough assumption that its value 
is the same as in the SM width. 

With this assumption, the relative difference between the production cross sections 
reported in [69] and a computation with the Higgs width on can reach the level of 2.5% 
in the high mass region, but is not significant for masses of mn < 160GeV, as shown in 
fig. 10. Let us remark, however, that in the high mass region the width of the Higgs boson 
is enhanced by the opening of new decay channels to third and fourth generation fermions, 
in a way that depends on the model. 

Ref. [69] refrained from providing cross-sections for Higgs boson mass values higher 
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Figure 10: Relative difference "~%wa • 100% of the cross section in SM4, within the two scenaria 
of [69], in the zero width approximation, a^^\, and in the approximation where T sMAij^h) — 
^SAiiTTiH)- Scenario 1 refers to a fourth generation down quark of mass m^^ = SOOGeV while 
scenario 2 assumes iridi — SOOGeV. In both cases the fourth generation up quark has a mass that 
is determined by eq.2 in section 3 of [69]. 

than TUfi > 300GeV. Cross-section predictions in that range should always take into 
account finite width effects. 



8. The Higgs cross-section for a variable bottom-quark Yukawa interaction 

In the standard model the bottom-quark Yukawa coupling is much smaller than the top- 
quark counterpartner, 

:^ = ^ ~ 0.02, 

and the bottom-quark fusion cross section is about 2 — 3% of the gluon fusion cross section. 
This feature may however not be conserved in extensions of the standard model, which for 
example contain more than one Higgs field electroweak doublets. In such a scenario more 
than one physical Higgs boson arise after electroweak symmetry breaking. The Yukawa 
couplings may then be modified by further mixing angles of the model and can differ 
strongly from their Standard Model values. It is then possible that the bottom-quark 
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fusion and gluon fusion via bottom-quarks become very significant in comparison to gluon 
fusion via top-quark loops. 

In this Section, we would like to study the sum of the two processes, gluon fusion and 
bottom-quark fusion, which contribute to the inclusive production of a Higgs boson as a 
function of the bottom Yukawa coupling (we denote by Y^ its value, normalized to the 
Standard Model). 
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Figure 11: Single Higgs production cross sections as a function of the re-scaling factor of the 
bottom-quark Yukawa coupling Yi,. The bands represent the uncertainty due to the choice of 
factorization and renormalization scale. 

In Figure 11 we demonstrate the inclusive cross-section for gluon fusion and bottom- 
quark fusion as well as their sum in the ZWA for a nominal Higgs mass of 120 GeV at the 
LHC for y^ = 7 TeV. We observe that for small Yj, the cross-section is dominated by gluon 
fusion. For high If,, bottom-quark fusion is dominant but there is also a large contribution 
from gluon fusion, however via bottom rather than top-quarks. Notice that the gluon 
fusion cross-section reduces for moderate values of Y^ due to a negative interference effect 
of top and bottom-quark loops. 

It is often the practice that uncertainties for beyond the Standard Model Higgs bosons 
are taken over from studies within the Standard Model. This may not be a very bad option 
if new physics only introduces new heavy particles and does not alter Higgs couplings to 
light quarks significantly, since all such scenaria can be well described by a common effective 
theory operator. In our scenario however, we need to be more attentive. Our calculation 
of the gluon fusion cross section is exact through NLO for both bottom and top-quark 
loop contributions. Our NNLO calculation includes only top-quark loops in the framework 
of HQET. For large Yj, where bottom-quark loops dominate, our evaluation of the gluon 
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fusion cross-section is reduced to NLO accuracy ( not NNLO). 




100 



Figure 12: Relative scale uncertainties of single Higgs production cross sections as a function of 
Vb- 

In Fig. 12 we demonstrate the scale variation of the cross-sections for gluon fusion , 
quark-bottom fusion and their sum, where we have combined scale uncertainties linearly. 
The scale uncertainty for gluon fusion increases for large lb as expected due to the dom- 
inance of the NLO only bottom-quark corrections. The scale uncertainty of the inclusive 
cross-section is dominated by the largest cross-section contribution. 

A large bottom-quark Yukawa coupling enhances the H ^ bb decay width. The 
total width of the Higgs boson can be derived from the Standard Model total width and 
branching ratios as: 



TniQ) 



-.SM 
H 



iQ) X [{Y,^ 



1 



BrT^^Q) + 1] 



.1) 



In Figure 13, we have plotted the combined bottom and gluon fusion cross section for 
a Higgs boson of mn = 120GeV for two different invariant mass windows. While the width 
effects are small up to Yf, ~ 30, they significantly change the total cross section for higher 
yb values, due to the steep increase of the Higgs width as the mass crosses the vector boson 
thresholds. 

9. The iHixs program 

The source code for iHixs can be downloaded from its website at 

http : //www . phys . ethz . ch/~pheno/ihixs 

Installation instructions can be found in the website and in the README file supplied in 
the distribution. Here, we briefly describe the main functionality of the code. 
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Figure 13: Combined gluon and bottom fusion cross section as a function of ; 
off-shell scheme. The bands represent the scale uncertainty. 



in on-shell and 



9.1 Usage 

The various features of iHixs are controlled by an input runcard, a text file that is edited 
by the user. To run with a given runcard as input type in the installation directory: 

./ihixs -i runcardjiame -o output_f ilename 

When no runcard is given, the program runs on the default card (called 'runcard') 
in the installation directory. When no output filename is given, the program writes the 
output in runcardjiame. out. 

The output consists of the total cross sections per perturbative order in QCD, together 
with the corresponding Monte-Carlo errors achieved and the PDF errors. Those are set 
to zero if no PDF uncertainty is requested in the runcard. The input runcard is also 
appended. 

9.2 Setting options and variables 

In the runcard anything after a hash symbol, '7^', is considered as a comment and is ignored. 
The following options are available: 

• pdf_provider : sets the PDF grid used. The user can choose between MSTW08, 
ABKM09 and GJR09. Within the MSTW PDFs there is also the option to switch 
confidence level from 68% to 90% and to use the MSTW grids with the strong coupling 
constant varied by one standard deviation from the best fit value. The exact filenames 
of the available grids are stated in the default runcard. 
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effective_theoryjElag : set to for the exact LO and NLO QCD effects and 
HQET approximation for NNLO. Set to 1 for the improved HQET approximation 
through LO, NLO, NNLO. 

no_error_f lag: Set to to calculate with PDF uncertainty, set to 1 to calculate 
without PDF uncertainty. 

collider: Set to 'LHC or 'TEVATRON' 

Etot: The total center of mass collider energy. This option is ignored if the collider 
chosen above is Tevatron. 

mhiggs: The nominal mass of the Higgs boson. 

higgs_width_scheme : Set to for the default finite width scheme. Set to 1 for the 
Seymour scheme, for a description of these schemes see section 4. 

higgs_width_grid: = The path^'' of the file with the grid for the width of the Higgs, 
and the branching ratios to 77, WW, ZZ and bb as a function of mn- If no path is set 
the default grid is used, HdecayGrid . dat, constructed with Hdecay v . 3 . 532 [26]with 
arguments that can be read in the header of the file. If the user supplies a grid file of 
his own, operating requirements are that the maximum number of grid points cannot 
exceed 16200, that the first three lines of the file are reserved for comments (so they 
are not read) and that the format of each line is respected, i.e. that the data is given 
in the order mn, ^h, BR^^, BRww-, BRzz, BR^j^. 

minjnh : Setting a minimum in the invariant mass of the Higgs boson. This allows 
the user to study the total cross section in the presence of kinematical cuts. 

maxjnh : Setting a maximum in the invariant mass of the Higgs boson. 

bin_f lag : Set to 1 to produce files with the bin-integrated Higgs invariant mass. 
Set to not to produce it. The data files produced contain the cross section per bin, 
with the bin size set to 1 GeV, from 30 to 2000 GeV at LO, NLO and NNLO. The 
files are named 'masshisto$r?7.jY.$order', so e.g. for mn = 200GeV the NLO file 
will be 'masshisto200. 1 '. 

muf /mhiggs : The ratio of the factorization scale and the Higgs mass. 

mur /mhiggs : The ratio of the renormalization scale and the Higgs mass. 

DecayMode: Set to no_width for the zero width approximation total cross section, 
to 'total' for finite width total cross section, or to the decay modes 'gamma gamma', 
'ZZ', 'WW, 'b b-bar'. 

ProductionMode: Set to 'gg' for gluon fusion or to 'bb' for bottom-quark fusion. 



^° Absolute or relative to the run directory. 
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• K_ewk: This is a global rescaling factor for all electroweak corrections. Set to 0.0 to 
switch them off. 

• K_ewk_real: Set to 0.0 to switch the electroweak corrections to H + j off. 

• K_ewkj:eal_b: Set to 0.0 to switch the electroweak corrections to H + j that include 
diagrams with massive quarks or Higgs boson in the loop, off. 

• m_top: the pole mass of the top-quark. 

• Gamma_top: the width of the top-quark. 

• Y_top : rescaling factor for the SM Yukawa coupling of the top. Note that this can 
be set to an arbitrarily small positive value, but not to 0.0 exactly. 



• m_bot: the MS mass of the bottom-quark at lOGeV. 

• Gamma_bot: the width of the bottom-quark. 

• Y_bot: rescaling factor for the SM Yukawa coupling of the bottom-quark. 

• heavy quark: Optional extra quarks in the model. The argument of this option 
should be formatted as niQ : Tq : Yq where Yq is the rescaling factor from a SM-like 
Yukawa coupling mq/v. For example, adding an extra 300GeV quark with width 
1.2GeV and a Yukawa coupling that is 5.7—^ the user should type: 

'heavy quark = 300.0 : 1.2 : 5.7' 

• m_Z: the mass of the Z boson. 

• Gamma_Z: the width of the Z boson 

• m_W: the mass of the W boson 

• Gamma_W: the width of the W boson 

• epsrel: Sets the relative Monte-Carlo integration error. 

• epsabs: Sets the absolute Monte-Carlo integration error. 

• nstart : Sets the number of points per Vegas iteration. 

• nincrease: Set the number of points by which the number of points per iteration 
increases 

• mineval: Set the minimum number of points before ending the Monte-Carlo integra- 
tion 

• maxeval: Set the maximum number of points after which the integration ends. 
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• adapt to central only: Set to to force Vegas to adapt to all integrand. Set to 1 
to adapt to the central integrand only. This is useful when running with PDF errors. 
Then each member of the PDF grid is treated as a separate integral. Adapting to 
the central only assumes that the peak structures of all integrals is similar which is 
a good approximation, and saves some CPU time. 

• vegas_verbose: Set to for silent Vegas output. Set to 2 to have information about 
each iteration printed in the standard output (the console). 

9.3 Libraries used 

The program uses the following libraries: 



• 



• 



• 



The Cuba library [78], v. 2. 1, for numerical integration. We use the Vegas algorithm, 
that employs importance sampling for variance reduction. For details on integration 
related arguments see [78] or the manual included in the Cuba-2 . 1 directory. We 
distribute Cuba-2 . 1 and compile it from source. 

The LHAPDF library [79]. We assume the library is installed by the user. See 
Ref [22] for details on installation. 

The package OneLOop [19, 20] for the evaluation of one-loop scalar integrals with 
complex masses. We use the library for the evaluation of finite box and triangle 
master integrals with massive propagators, necessary for the electroweak corrections 
to H + j with massive fermions in the loop. We have checked our implementation of 
all other integrals against both OneLOop and QCDloop [21] at the limit of zero width 
for the massive propagators. We distribute OneLOop and compile it from source. 

The CHAPLIN package [18], for evaluating harmonic polylogarithms up to weight four 
for any complex argument. This package is also distributed. 



10. Conclusions 

In this article, we have presented a computer program, iHixs, for the inclusive cross-section 
of the Higgs boson in gluon fusion and bottom-quark fusion. iHixs provides the most 
precise predictions for the Higgs boson rate at hadron colliders in fixed order perturbation 
theory, including QCD corrections through NNLO and electroweak corrections for virtual 
and real radiative partonic processes. iHixs is interfaced with the LHAPDF library and 
allows the assessment of uncertainties due to the various determinations of the parton 
densities which are available in the library. 

In a time in which Higgs boson searches are growing in intensity, iHixs provides a very 
flexible tool which can assist in this effort. iHixs provides inclusive cross-section predictions 
in and beyond the Standard Model, by allowing modifications of Yukawa and electroweak 
couplings as well as the introduction of new quarks with arbitrary Yukawa couplings and 
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masses. In addition, one can readily introduce effective Higgs-gluon interactions whicli can 
account for further beyond the Standard Model effects ^^. 

The phenomenology of the Higgs boson and its production rates have been described 
extensively in the literature (recent updates can be found in Refs [50, 73-75]). We ded- 
icated a very short analysis to issues which have been studied earlier at length, such as 
the magnitude of perturbative corrections and the convergence of the perturbative series 
beyond NNLO in QCD [76]. 

We have noticed that in recent experimental studies [3] the zero width approximation 
which is used for an expected light Higgs boson in the Standard Model is also employed for 
Higgs bosons or Higgs bosons with a sizable decay width. In this publication, we discuss 
finite width effects on the cross-section due to resonant Higgs boson diagrams. We also 
employ a prescription to estimate the effect of the signal-background interference for high 
Higgs boson masses which can be dramatic. 

We believe that a realistic description of the Higgs line shape is necessary in setting 
exclusion limits for the Higgs boson. We remark that the description of the line shape in 
parton shower Monte-Carlo generators can be very different (for a comparison see Ref. [77]). 
We have also demonstrated that the magnitude of radiative corrections (K-factor) differs 
from expectations in the zero width approximation when a large range of virtualities for 
the Higgs boson is sampled. 

iHixs allows the user to perform exhaustive studies of the Higgs boson cross-section at 
hadron colliders. We are looking forward to comparing iHixs predictions with LHC data. 
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A. Tables of Higgs cross-sections 
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-3.73 
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5.89 
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Table 5: Total cross section for LHC at y/s = 7TcV with MSTW PDF errors (corresponding to 

68%CL). 
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4.1 


-4.1 


6.13 


-8.57 


260.0 


2.77 


4.2 


-4.2 


6.03 


-8.52 


270.0 


2.58 


4.3 


-4.3 


5.94 


-8.45 


280.0 


2.41 


4.4 


-4.4 


5.82 


-8.4 


290.0 


2.27 


4.5 


-4.5 


5.7 


-8.35 


300.0 


2.15 


4.6 


-4.6 


5.55 


-8.28 



Table 6: Total cross section for LHC at ^/s — 7TeV with ABKM PDF errors (corresponding to 

68%CL). 
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niH 


a{pb) 


%Spr)P 


%Sppp 


%v. 


%^u 


110.0 


18.66 


3.6 


-3.6 


7.87 


-8.63 


115.0 


17.1 


3.5 


-3.5 


7.73 


-8.59 


120.0 


15.79 


3.5 


-3.5 


7.58 


-8.55 


125.0 


14.58 


3.5 


-3.5 


7.46 


-8.5 


130.0 


13.49 


3.4 


-3.4 


7.35 


-8.46 


135.0 


12.52 


3.4 


-3.4 


7.25 


-8.42 


140.0 


11.66 


3.4 


-3.4 


7.14 


-8.39 


145.0 


10.88 


3.4 


-3.4 


7.04 


-8.36 


150.0 


10.17 


3.3 


-3.3 


6.95 


-8.33 


155.0 


9.53 


3.3 


-3.3 


6.85 


-8.31 


160.0 


8.89 


3.3 


-3.3 


6.77 


-8.3 


165.0 


8.09 


3.4 


-3.4 


6.73 


-8.29 


170.0 


7.53 


3.4 


-3.4 


6.66 


-8.31 


175.0 


7.05 


3.4 


-3.4 


6.58 


-8.33 


180.0 


6.62 


3.4 


-3.4 


6.53 


-8.36 


185.0 


6.17 


3.4 


-3.4 


6.46 


-8.42 


190.0 


5.78 


3.5 


-3.5 


6.4 


-8.46 


195.0 


5.45 


3.5 


-3.5 


6.35 


-8.5 


200.0 


5.16 


3.6 


-3.6 


6.29 


-8.54 


210.0 


4.66 


3.6 


-3.6 


6.18 


-8.62 


220.0 


4.24 


3.7 


-3.7 


6.04 


-8.7 


230.0 


3.88 


3.8 


-3.8 


5.94 


-8.75 


240.0 


3.57 


4.0 


-4.0 


5.83 


-8.81 


250.0 


3.31 


4.1 


-4.1 


5.7 


-8.9 


260.0 


3.08 


4.2 


-4.2 


5.61 


-8.92 


270.0 


2.88 


4.4 


-4.4 


5.49 


-8.96 


280.0 


2.72 


4.5 


-4.5 


5.38 


-8.97 


290.0 


2.58 


4.6 


-4.6 


5.24 


-9.0 


300.0 


2.46 


4.8 


-4.8 


5.09 


-9.01 



Table 7: Total cross section for LHC at -^/i 
68%CL). 



7TeV with GJR PDF errors (corresponding to 
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B. Matrix element coefRcients 

B.l Master Integral Definitions 

In the following we will use the shorthand notation 

Pi\i2--in ^ Ph ' Pi2 ' •• 'Pin 

and define the Mandelstani variables 

s = {Puf, t = (P23)^ u = (P13)^ rriH = {Pl23f 
where Pi = P2 = P3 = 0. The master integrals are then given by 

d^k 1 



Tadp(m 



2\ 



i7r^/2 [fc2 - m2] 



iiuh[s,m 



Tria(s,m ) = s ■ 



i7r^/2 [/i;2 _ j7j,2][(^/j _^ p-^^'^2 _ j^2 

d^k 1 



ivr-D/^ [A;2 _ j^2j 1^^^ + pi)2 — m2] [(fc -(- ^12)^ — w-^] 

2n /" ^^^ 1 

Box(s, t,u,m ) = s ■ t ■ ■ 



i7r-D/2 [fc2 — 77i2] [(/;; + pi)2 — m'^] [{k + pi2)2 — m2] [(/j + pi23)2 — m^] 



-1—1/ 9 9\ I Qi rv X 

BubE(s,m]^, 771-2 ) 



77r^/2[fc2-m2][(A:+pi2)2-m2] 



TriaE(s,?77|^,mf,m^) = s- / —^tt^ 2T7m ^2 STTTT- T2 2? (^-1) 

77r^/"' [A;^ - mf J [(/c + P12) - m^\ [{k + P123) - ""^iJ 



TriaF(s, 777]^, 7/1,2) ~ ^' 



d^k 1 



1—1 /999\ I Qj rC L 

BoxE{s,t,mH,m„m,) = s-t- J ,^d/2 [p _ ml][ik + p,y - ml][ik + pu)' " ml][ik + p^23? - m?] ' 

We also define 

BubD (s, m ,tj = Bub (s, 777„) — Bub (t, 777 ) . 



B.2 A 



gggH 

2^ 



, 2 (2 St + t77 + STt) Box (s,t, 77, 7772) (2 t77 + St + STx) BoX (t, 77, S, 
KaaH (S' *' ^' "ij = 2 — + 2 : 



77 S 

(2 SU + tU + St) Box (77, S, t, 7n.2) (4 t77 + ^2 + 7^2) sBubD (s, 7n,2, S + t + 7i) 

(7*2 + 4 S7t + s2) tBubD (t, 77l2, s + t + 7i) (4 St + S^ + t^) TiBubD (7i, 777,2, s + t + 77 



\2 



+ - 



(77 +S)^ {t + s) 

5 S'^7it2 + 5 S7i4^2 _,_ g g^S^S _^ ]^q 52^8^2 _^ jO 527(2^3 _^ IQ s37i2t2 + g 587^^3 _^ 5 ^4^2^ 

t (77 + t) S (77 + s) 77 (t + s) 
5 sH'^U + 5 ^527^4 _^ g ^^3^3 _^ 5 gi^2^ j^ 2 77^t3 _^ 3 77^^"^ + 2 5^77^ + 2 s'^Ti^ + 2 sH'^ + 2 s'^t^ 



+ 4 



t{u + t) s{u + s)u{t + s) 
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-16m„ 



j.3„,3 I Qj-2„,4 I Qj.4„2 



ru^s + 4 tu^'s + 8 tn^s^ + 8 s^f^w + 8 t^u^s + 4 Tsu + 4 f^w^ + 3 t^n* + 3 rs 

{u + sf {t + s)'^ {u + tf 



+ - 



3 s^u'^ + 3 t*^u2 + 4 s^u^ + 3 s^'n^ + 3 sH^ + 4 t^s^ + 8 s^uH + 4 s^nt + 6 sH^u^ + 8 s^t^n' 

Tria (s + t + u, m^) 

u^ + t^) ml 2 t^s + 3 st^u + 3 suH + 2su^ + t^u + u^t 



+ 16 



16 



+ 16 



(u + t) 



tu {u + t) 



[s^ + u^) m^ 2 ts^ + 3 s^ut + 3 su^t + 2uh + su^ + s^u 



iu + s) 



4.2 , „2\ ^2 



SU [U + s) 



[t^ + s^) m^ 2s^u + 3 s^ut + 3 st^u + 2t^u + t^s + ts^ 



it + s) 



st (t + s) 



I Tria (s, m^) 
Tria (t, m ) 
Tria (n, m ") 



-16 



(n + t){u + s) {t + s) 



^llgHi^^'t^u.ml) = [-lQml-2— ^ - 

2 st + s"^ — 2tu 



Box (^s,t,u,mg) 



+ ( -16m^ + 2- 



Box (t, u, s, m^) 



+ 16m2+4— )Box(u,s,t,m2) +8— \^.^1 1 



+8 



(m^ + 4 sti + s^) tBubD [t, nig, s +t + u) {s - t) uBubD {u,mg,s + t + u) 



+ ( -16 
-4 



{u + sY t + s 

[uH + 3 t^u^ + 5 s^ut + 3 sH"^ - su^ + b suH + 8st^u- s^u + ts^) m^ 



{u + sY {t + s){u + t) 

2 sH"^ - st^U^ + 2 t^t^S _ ^2^2^ _ 2 i^2^2 _ 2 s2^3 _ 2 ^3^2 



+ 16- 



(it + s) stu 
-u + t) ml ,2t^s-t^u-2u^s 



Tria (s + t + u, m^) 



li + t 



+ 4 



tu 



Tria (s, m^) 



+ 16 



{u^ + A SU + s^) m\ ,2ts^ -2s'^u^ + s'^ut + su^t + 2uh 



+ 4- 



Ui + s 



SU {u + s) 



Tria (t, m^) 



+ -16 



(s-t)m^ 2s'^u + t^s-2t^ 



t + s st 

^llgHi^'^^'^^'^l) =A2it,U,S,ml) 
^33gi^(^'*'^'"^g) = A2{u,S,t,ml) 
B.3 AqqgH 

Defining 



Tria (u, m,^ — 16 



n. SU+ St + tu 

^ \ 1 ^ 



u + s 



(B.2) 

(B.3) 
(B.4) 



4m^ 
^1 = 32' 



m 



H 
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we have 

3 / -2A{Tq 



^qqgH {Tq,y,m]j) 



4] 1-y 

Tvia{y ■ m^,m'i^ 



+ (1_^^)2 ^^^P(^ • "^H, rnl,ml) + j^^^ 



y 



{l-e)Tx\a{y-rr?H,ml)) (B.5) 



^_y.- -.—-.. .-.,■■, ^ 
Notice that for yi2 — )■ we have that 



Wm. Aqqg {Tq,yi2,mjj) = A {Tg,mjj) , (B.6) 

J/12— )-0 

which is the familiar kernel of the born gg ^ h amplitude. 

B.4 Aewk 

f ( TYl S 1 

Aewk (s, t, u, mH,mz) = r 2( \+\_ j-1 BoxE("'='^^ (s, t, m\, ml,0) 



m 



(m^ - s) 



"^ — i — 77 1^ — ^TriaE (s, m^;^, m 0) 



1 

+ - 

s 



1 



2 4 

mt rnt 



s + u mi is + t) — st 



TriaE (t, rriu^m^, O) 



, Bub(m^,m^)-BubE(t,m^,0) 
s{s + u) 

with u = rnjj — s—t. Where we have expressed the form factor in terms of the 6-dimensional 
Box using the following relation 

BoxE(-^) {s,t,ml,ml0) = - ^^' + '^ "^^^ 'f^ ^^' \') ^f " f} ^"^ + f ^^^^ ^^^'^ 

{{s + t) m^^ - st) BubE (s, m^^ O) ((s + t) m^^ - st) BubE (t, m^^ o) 

2tu {rriz^ ~ s) e 2su {niz^ ~ t) ^ 

[[st + t'^ - SU + tu) rriz'^ - st{u + t)) TriaE [s, niH'^^mz'^, O) 

2{-l + 2e)stu 
[[su + s'^ - tu + st) ruz'^ — st{s + u)) TriaE [t, m/j^, m^^, O) 



2{-l + 2e)stu 
[s + tf m/ -2st{s + t) ruz"^ + sHA BoxE (s, t, mn'^^mz'^, O) 



2(-l + 2e)stn ^ 



B.5 AZ, 



The form factor for the mixed QCD-electroweak corrections to H + j with a top-quark in 
the loop is: 

, 2 2\ "T-t^Tadp [niup') - mt^Tadp (mj^) 
AZk{s,t,u,ml,mi) = ^ ^^^^4 ^ 
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+mt^ i-mui^ + mt^) BubE (s, niui^^mt^) 



+ 



+ 



[—mu'^mt^ — 2 rrit^mw^ — 4 m^^) BubE (t, m^^, mt^) 

4s (u + s) niw'^ 
[mn'^mt^ + 2mt^mw^ + 4mi„^) Bub {niH'^ ,mw^) 



As {u + s) rriu,'^ 
+ \ [- 2t{-t + s)mw^mt^ -t{-t + s)mH'^mt^ - 2t{-t + s)mj^mt'^ 

+t (3 s + us — tu — t ) niyj mt —t smn rnt + At {—t + s) rriw mt 

-2 1 (s^ - st) rut'^myj^ + 2 s^t^m^^mj^]TriaF (t, niw^ ,mt^) 

+ [- {s + t) mn'^mt^ + (-2 1 - 2 s) mj'mt^ - m//^ (tn + us + t^) mt^ 

+ (t^ + 3 s^ + 6 st + tu + ns) mj'mt^ 

+ (2 s + 2 1) m^'^mt*^ + mi:^ ^ (stu + t^s) rrit'^ 

+ (tu^ + t^ + 3 s^u + 3 stn + 2 t^s + 2 t^u - s'^t + su^) mj'mt'^ 

+ (—5 1^ — 7 s^ — 5 ns — 5 tti — 8 st) mw'^mt'^ + (6 1 + 6 s) m^^mt'^ 

—2s t{u + t) niyj mt + (—3 s t — 3 stn — 2s n — 3t s) m^t, ttij 

+ (l2st + 5tu + 5us + 9s^ + 5t^)m^'^mt^ + (-10 s- lOt) m^^mj^ + 4m^„^s^t 

+ {-^st-As^)mJ' + (4t + 4s)m^^°]TriaE(s,mH^,r7T,^^,mt^) 

+ [mt^{s + tf + (-s^n + 3t^s + 2s^t + 5stu- s^) mt^ 

+ (-4 s^ - 2 1^ - 2 st) m,Jmt^ - s^t (n + s) m^^ 

+ {s^u + ?,s^ -stu-2sH-?> t^s) mjmt* + (-s^ + t^ + 4 st) my^^nit* 

+2 s^tmyj'^mt^ - 2 s^ {-t + s) m^^'mt^ + (-4 st + 4 s^) m^'^mt^]TriaF (s, mw^,mt^) 

+ [2 1 (s + t) mjmt^ +t{s + t) rxiu^rat^ 

—2t(^st + tu + t +s — us) mt m,yj +tmH {st + 2us)mt 

-6 1 (s + t) mj'mt^ + t (5 s^ + t^ + tu + 9 ns + 8 st) mt'^m^,'^ 

-t [t'^s + 4 s^u + s^t + stn) mt^myj"^ + At {s + t) mj^ - 4 1 (2 st + s^) 

+4m^„^s2t2]BoxE('^=^) (s, t, mn'^ ,mj ,mt^) 



m,„^ 



X Asm^'^ yAmt^stu + ((s + t) {-m.J' + mt^) + st) 

+ <^ (s + t) (2 s + n + t) mh^mt^ + 2 (s + t) (2 s + u + t) m^Jmt^ 

-2 {s + t){2s + u + t) mw*mt^ 

+ (-3stu - gs^t - 4s^ - su^ - Qt^s - 3t^ - 4t2n - tu^ - Ss^n) mjmt^ 

+ (t^n + 3 t^s + 6 stn + tn^ + 3 s^t + sn^ + s^ + 2 s^u) m^'^mt'^ 

-6 (s + t) (2 s + u + t) m.Jmt^ 

+ (17 s\ + 5 sn^ + 27 stn + 5 tn^ + 6 1\ + 17 s^t + 12 t^s + 12 s^ + 3 1^) m^^'mt^ 

+ (-14 t^sn - 6u^st -2t^u- 9 sh^ - 3t^u^ - tu^ - 7 s^u -5s^u^ - 6t^s 
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-su^ - 8 s^t - 3 / - 13 s'^ut) mjmt'^ + (2 sH"^ + t^su + u^st + 2 s^nt + sH) mn^mt^ 

+10 (s + t) (2 s + n + t) mj'mt'^ 

+ (-19 s\ - 5 su^ - 41 stu -5tu^ -8 t^u - 33 s^t - 22 t^s - 14 s^ - t^) raj'mt^ 

+ (t^n + t\^ + 12 t^su + 15 s^i^ + 12 n^st + 2 s^v? + 3 t^s + 4 s^n 

+ 10 s^t + 2 s^ + 22 s^ut) ra^^nit^ 

-st (2 sn^ + 2 s^ + t^u + tu^ + 4 s^n + 3 stu + 2 t^s + 2 s^t) m,J'mt^ 

-A{s + t){2s + u + t) mj^ + (4 t^n + 12 stu + 4 s\ + 20 s^t + 12 t^s + 4 s^) mj 



2 2\ 



-4 stm^^ (2 us + 3 st + 2 s^ + 2 tn) + 4 s^t^ (n + s) m^„^ llViaE (t, mu^ , 

X 4s(u + s)m^^ C4mt^stn+ ((s + t) (-m^^ + mt^) +st)^^ (B.9) 

where 
BoxE('^=^) (s, t, m/f ^ m^^ mt^) = ish'^ - 2 (s + tf mt^m,J - 2 st (s + t) m.J' 
+2 st (t + s + 2 u) mt^ + (s + t)^ m^^ + (s + t)^ rrit 



2 2 4 , / , ,^2 4^ BoxE (s, t, m/f 2, m^^^ ^^2^ 



+ 



+ 



2stu(-l + 2e) 
[[su — tu — st — t^) mw^ + (tu + st — su + t^) mj^ + st (u + t)) TriaE (s, m//^, m^^, m^^) 

2stu(-l + 2e) 
{{—st — su + tu — s^) m^^ + {su + st — tu + s^) ?n,t^ + st (s + u)) TriaE (t, m//^, m^^, m^^) 

2stu(-l + 2e) 
(sTriaF (s,m^^,mj^) + tTriaF {t,mw^ ,mt^)) ((-t — s) m^^ + (s + t)mt^ + st) 



2stM(-l + 2e) 



(B.IO) 



and m?^ = s + t + u. The ttij — )■ limit of this form factor trivially lead to Ag^fc of 
eq. B.7. 
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